Poisson transition rates from time-domain measurements with finite bandwidth* 
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In time-domain measurements of a Poisson two-level system, the observed transition rates are 
always smaller than those of the actual system, a general consequence of finite measurement band- 
width in an experiment. This underestimation of the rates is significant even when the measurement 
and detection apparatus is ten times faster than the process under study. We derive here a quan- 
titative form for this correction using a straightforward state-transition model that includes the 
detection apparatus, and provide a method for determining a system's actual transition rates from 
bandwidth-limited measurements. We support our results with computer simulations and exper- 
imental data from time-domain measurements of quasiparticle tunneling in a single-Cooper-pair 
transistor. 

PACS numbers: 02.50.Ga, 85.35. Gv, 72.70.+m, 85.25.Cp 
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We consider here a physical system switching inco- 
herently between the states A and B in an alternating 
Poisson process, with characteristic rates V a = t7 and 
Tb = t~b i where ta and tb are the lifetimes of the re- 
spective states [Fig.^a)]. This two-state model has been 
used in the analysis of a wide variety of problems arising 
in medicine reliability theory 0, network traffic 0, 
cell physiology 0, 0| , materials science 0, , and con- 
densed matter physics H 0] , especially in relation to 
1// noise Time-domain measurements of such sys- 
tems produce random 'telegraph signals' that represent 
the underlying transitions in the system convolved with 
the response of the measurement apparatus [Fig. ^b)]. 
The transition rates between the states Ta and T b often 
contain information about the underlying physical mech- 
anism 0, El E3 , and one of the goals of the experi- 
ment is to extract these rates from the data. To do so, 
one employs a detection algorithm that operates on the 
measured telegraph signal |3, E3- El and, by means 
of a statistical test, determines the dwell times in each of 
the states between transitions, Fig.^c). The dwell times 
are then histogrammed to give the lifetime distribution 
in the two states, from which the transition rates may be 
determined. 

A comparison of Figs, da) and (c) reveals that what 
is observed in an experiment, i.e., the state of the de- 
tector, does not always reflect the true state of the sys- 
tem. Particularly, the bandwidth of the measurement 
limits the visibility of the underlying process on short 
time scales. In this Letter we show that the response 
of the measurement chain modifies the statistics of the 
observed process, and that the experimentally obtained 
transition rates always underestimate those in the under- 
lying system. We argue that this is a general feature in 
finite-bandwidth time-domain measurements of stochas- 
tic two-state systems, which introduces systematic errors 
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FIG. 1: Simulated data shown at different stages of the mea- 
surement process: (a) The true state of the system, (b) the 
observed telegraph signal, and (c) the detected state sequence. 



in estimating the actual transition rates from the data. 
Using a straightforward model, we calculate the lifetime 
distribution in the observed process, relate the experi- 
mentally observed transition rates to those in the under- 
lying process, and compare the results of our analysis 
with both experimental and simulated data. 

In the following we assume that the overall response of 
the measurement chain and the detection algorithm may 
be characterized by an effective detection rate, Tdot = 
1/Tdetj where Tdet is the mean time taken by the detector 
to register a transition in the measured observable after 
a transition in the underlying system has occurred |l5| . 
It depends on the physical bandwidth of the signal, the 
details of the detection algorithm, and the signal-to-noise 
ratio of the measurement. The statistical nature of a 
detector operating on a stochastic (and usually noisy) 
signal allows us to treat the detection itself as a stochastic 
process. To simplify the discussion, we will assume here 
that the detection process is homogeneous, so that Tdet 
is independent of time 0] . 

We note that the detection problem has received con- 
siderable attention in the context of patch-clamp record- 
ings on cellular ion channels [Illl^ll&llJI- These authors 
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FIG. 2: Schematic state diagram for the detection process. 

have tended to characterize the detector with a constant 
delay time, an assumption that we argue is not physical 
and does not agree with our experimental results shown 
below. Their probabilistic approach to the problem, how- 
ever, allows for generalizing the analysis to processes with 
extended state spaces. 

To understand how the finite bandwidth of the mea- 
surement modifies the statistics of the observed process, 
let the underlying process be Poissonian, so that tA (is), 
the dwell time in the state A (B), is an exponentially 
distributed random variable with mean ta (tb)- In a 
particular sequence, A — > B — > A, the excursion into B 
might not be registered by the detector if ts < T dct ■ Two 
consequences follow: 1) the lifetime distribution will have 
a short-time cutoff near Tdct! and, more seriously, 2) the 
two neighboring occurrences of A, which were originally 
distinct, will be effectively 'stitched' together to register a 
single event, longer in duration. Therefore, the observed 
distribution of t& will be artificially reweighted towards 
longer times, with significant consequences for inferring 
the system parameters from the data. We will calculate 
this new distribution from the dynamics of the detection 
process, using the model shown in Fig. [21 

In our model we separate the occurences of the states 
A and B into two manifolds, which we condition upon 
the state of the detector - either A* or B*, according 
to the state into which the last detected transition has 
occurred. Suppose that the detector has just entered A*. 
We are thus on the left hand side of Fig. [3in the state A. 
The system may then switch into state B, from which a 
transition into the B* manifold (detection) will occur at a 
rate T^ct- While waiting to be detected, the system may 
cycle several times between A and B until a transition is 
eventually registered. Once in B* the cycle starts again 
until we return to A*, and so on. 

The observed process is thus A* ^ B* . To find the 
distribution of times in A* before a transition is detected 
(A* — > B*), we write the rate equations that describe 
the evolution of the system in the left hand side of Fig. [21 
(the dwell-time distribution in B* may be found by in- 
terchanging A and B in the following), 

a* f Pa = -TaPa+T b Pb, m 
\Pb = -{Tb + T dct )P B + T A P A . [ ' 



The experimentally accessible quantity, the dwell-time 
histogram, is the probability density of leaving A* in 
the interval [t,t + dt], which is proportional to h(t) = 
rdet-Ps(i)- Solving Eq. Q with the initial conditions 
Pa(0) = 1 and P B (0) = 0we obtain 

h(t) = | T A T dct e-^ xt wah(0t/2), (2) 

where \ = T A +T B + T deU and 6 = VA 2 - 4T_4r dct . A 
similar expression has been obtained in Ref. 19]. 

Eq. (J2J is our main result, and has several important 
implications: 1) the experimentally observed process is 
no longer Poissonian. This should be kept in mind, for 
example, when the correlations or the full counting 
statistics |2lJ of the observed process are analyzed. In the 
short time limit, where h(t) approaches zero, the ability 
to observe the state B is bottlenecked by Tdct- In the 
long time limit, however, Eq. (J2J) tends to an exponential 
distribution. 2) The mean dwell time in A* is (tA*) — 
J th(t)dt — A/r^Fdct, while the rate parameter of the 
long-time exponential tail of h(t), which we call Ta*, is 
given by 

t a - = \ (i - - 4r A r dot /A 2 ) . (3) 

3) Researchers have traditionally interpreted both exper- 
imentally derived quantities, the lifetime ta* = 1/Ta* 
and the mean dwell time (tA*) (which are not equal), 
as the Poisson lifetime ta, but neither quantity is an 
accurate estimate of the true lifetime of the state A. 
For example, even when the measurement bandwidth 
is an order of magnitude greater than the process, say 
ta = tb = 10id ot , we see that ta* and (tA*) overes- 
timate ta by 10 % and 20 % respectively. In the limit 
Tdet/TA - > 0, the fractional correction to the lifetime 
vanishes as T dct /TA- 4) Lastly, if (but only if) the de- 
tection rates are equal for both states, then Eq. (J2J) 
also implies that the equilibrium occupation probabili- 
ties in the observed and underlying processes are equal, 
Pa* — (tA*) / {(tA*) + (tB*)) = Pa- This somewhat sur- 
prising result also follows from the steady-state solution 
of the master equation corresponding to Fig. [21 

We apply the results of the above analysis, mainly 
Eq. J2J, to an experiment 22] in which we measured 
the dynamics of quasiparticle (QP) tunneling in a single- 
Cooper-pair transistor fSCPTl El E^ . The tunnel- 
ing of quasiparticles onto and off of the island of an SCPT 
is expected to follow a Poisson process, taking place on 
microsecond time scales. The transistor was operated 
in a regime where the island can trap a single QP [25| . 
so that and Tb correspond respectively to the QP 
emission (un-trapping) and capture (trapping) rates. To 
observe this process, we used an rf rcflcctometry tech- 
nique with the transistor held at zero dc bias |27j . 
The reflected rf signal, whose magnitude indicates the 
presence or absence of a QP on the island, was measured 
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FIG. 3: (a) Histograms of QP dwell times on the island 
of the SCPT. The three curves correspond to different IF 
bandwidths- 0.3, 1, and 3 MHz, and are offset for clarity, (b) 
Measured lifetimes vs. effective detection time for state A (QP 
present, main panel), and state B (QP absent, inset). Squares 
and triangle represent analysis using the Schmitt-trigger and 
Page-Hinkley algorithms, respectively. The solid lines were 
calculated using Eq. and its equivalent for Ts*, with 
ta = 26 /is, tb = 3 fjs, and circles show the results of our 
simulations (see text). 



with a spectrum analyzer. The video output of the ana- 
lyzer was digitized at 500 ns intervals and recorded by a 
computer, where a typical record followed the evolution 
of the system over a span of one second. 

For this analysis, we recorded a series of time traces 
that differed only in the bandwidth of the intermediate- 
frequency (IF) filters of the spectrum analyzer. Each of 
these time traces was analyzed by use of both a simple 
hysteretic ('Schmitt trigger') detector algorithm 0], and 
a Page-Hinkley cumulative likelihood ratio algorithm 0, 
ll^ . In both cases we set the detection thresholds so 
that the detector's false-alarm rate was negligible, at the 
expense of compromising the detection efficiency. 

In Figure |2Ia) we show histograms of QP dwell-times 
on the island of the SCPT. The histograms were ob- 
tained using the Schmitt-trigger algorithm to process 
time traces measured with IF bandwidths of 300 kHz, 
1 MHz, and 3 MHz. It is clear from the slopes of the 
measured distributions (on a log scale in the figure) that 
the observed QP tunneling rate depends on the measure- 
ment bandwidth. It appears greater when measured at 
a higher bandwidth, in accordance with our understand- 
ing of the detection process. Also note the sharp drop in 
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FIG. 4: Observed lifetimes ta* from a simulated process 
with ta = 25 /xs and tb = 1 to 25 fis. The signals were 
filtered with f c = 1 MHz and processed with the Schmitt- 
trigger algorithm. From a fit (solid line) we find Tdet = 1.8 /xs. 
Because Tdet represents the overall response of both the filter 
and the algorithm, it is not surprising to find Tdet > !//<=■ 



the number of registered counts at short times, reflect- 
ing the detection bottleneck evident from Eq. J2J). This 
short-time behavior of the dwell-time distribution might 
be overlooked if the data are binned too coarsely. 

We obtain the QP tunneling rates, as measured with 
each bandwidth-algorithm combination, by fitting the 
lifetime distributions to an exponential. The experimen- 
tal values of the lifetimes ta* = 1/Ta' (QP on the island) 
and tb* = 1/T.B* (QP off the island) are shown respec- 
tively in the main panel of Fig. IHIb) and in the inset, 
plotted versus the effective detection time of the corre- 
sponding measurement. To estimate the detection time 
Tdet directly from the data, we chose to use the time 
at which the lifetime distribution peaks. This is only 
an approximation, giving an uncertainty in T^ ct of about 

±1 yliS. 

The true lifetimes are estimated by extrapolating the 
measured lifetimes in Fig. |3^b) to T^ct — 0, and are ap- 
proximately ta — 26.0 ± 0.7 /is and t b = 3.0 ± 0.5 /is. 
The solid curves in the figure were calculated from Eq. @ 
and its equivalent for Tb*, using the above numbers for 
t~a,b- We have also performed computer simulations, 
whose results are shown in Fig.^b) as circles. We gen- 
erated a series of ~ 10 4 transitions between two signal 
levels, whose durations were taken as random variables 
from two exponential distributions with means equal to 
ta.b above. We added Gaussian noise to the signal, and 
applied to it a range of convolution filters with varying 
bandwidths. The resulting signals were then processed 
using our Schmitt-trigger algorithm, and analyzed as de- 
scribed in the preceeding paragraphs. The agreement 
between the simulations, experiment, and the theory, ev- 
ident from the figure, indicates that our model captures 
the main features of finite bandwidth time-domain mea- 
surements. Characterization of the detector with a con- 
stant "dead-time" as in Refs. H3, E3| 

would result in an 

exponential dependence of the observed lifetimes on the 
bandwidth of the experiment, in disagreement with our 
results. 

Another set of simulations is shown in Fig. Here we 
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programmed a waveform generator to give an analog sig- 
nal that follows a simulated Poisson process. This analog 
signal was used to amplitude-modulate the rf carrier in 
our experimental setup. We varied tb while holding T4 
fixed, and used a fixed filter. We plot in Fig. 0] the ob- 
served ta* as circles, and the solid line is a fit to Eq. © 
with Tdct the only free parameter. A similar procedure 
may be used to determine the detection rate in a given 
experiment. 

Since only the observed process is accessible in an ex- 
periment, we proceed to find a transformation that ex- 
presses the true rates of the system in terms of the ob- 
served ones by inverting Eq. © and the correspond- 
ing equation for Tb* ■ The inverse transformation be- 
low is given in terms of the dimensionless quantities u = 

r^/rdetj v = r B /r dct , u* = rv/r dct , v* = r B ./r dct , 

and is easily generalized to the case of asymmetric detec- 
tion rates, 



1 — u* — v* 
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By use of Eq. J2J it is straightforward to analyze the 
propagation of uncertainties from the measured lifetimes 
to the true-lifetime estimates. The lifetimes in the un- 
derlying system can alternatively be extracted from the 
moments of the dwell-time histograms, Eq. : ta = 
T dct ((*>!*) 2 ~ 2 (*A*))> with a similar equation for tb- 

Although the example analyzed in this Letter is from 
the realm of condensed matter physics, we stress that our 
approach to the detection problem is applicable also in 
other contexts. There are a number of examples in the lit- 
erature for the treatment of imperfect detectability of an 
underlying Markov process in medical studies [28( , and of 
Markov processes that are only partially observed, e.g., 
in reliability theory and criminology |29l | . These authors, 
however, considered the detection problem as static and 
independent, by assigning the detector a fixed fidelity. 
Accounting for the dynamics of the detection process, as 
we have done here to address the effects of a detector's 
finite response time, may help identify and correct sys- 
tematic errors arising in studies in those fields. A more 
accurate description of the detection process has been 
developed in the field of cellular physiology 
however, our results suggest that our stochastic model 
of the detector is more realistic than the "constant dead 
time" model commonly used in that context. 

To conclude, we have shown that when a two-state al- 
ternating Poissonian system is measured in the time do- 
main, the statistics of the experimentally observed pro- 
cess depend on the response time of the detector. Partic- 
ularly, the observed process is no longer Poisson, and the 
transition rates that one obtains from the data always 
underestimate those in the underlying process. We have 



given analytic expressions for the relation between the ex- 
perimentally observed rates and their actual values, and 
shown that our results are in good agreement with both 
experimental and simulated data. We argue that not ac- 
counting for the effects of finite measurement bandwidth 
will lead to results that are inaccurate at best, even when 
the bandwidth of the experiment is an order of magni- 
tude greater than the underlying process. On the other 
hand, measurements of fast processes approaching a well 
characterized bandwidth limit of an experiment can still 
give meaningful results using the analysis presented here. 

We thank R. L. Kautz, S. Nam, and J. A. Stroscio for 
valuable discussions. 
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